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DESCRIPTION 

RADIATION IMAGING DEVICE FOR CORRECTING FOR BODY 
MOVEMENT, IMAGE PROCESSING METHOD, AND COMPUTER 
PROGRAM 

TECHNICAL FIELD 

The present invention relates to a radiation imaging device, an image 
processing method for the radiation imaging device, and a computer program 
for causing the radiation imaging device to execute the image processing 
method. In particular, the radiation imaging device and the image processing 
method according to the present invention can be suitably used to reduce in 
radiography an artifact which occurs in a tomographic image due to movement 
of a body to be tested (this movement is called "body movement", hereinafter). 

BACKGROUND ART 

Fig. 9 is a diagram showing a concept of a cone beam X-ray CT 
(computed tomography) device in which a plane sensor is used. 

In Fig. 9, » X-rays * N , s h! -^w from an X-ray source 901 is 
transmitted through a body 903 to be tested, and the transmitted X-ray s are i-s- 
detected by a plane sensor 902 provided at the position opposite to the X-ray 
source 901. Here, it should be noted that the body 903 to be tested is also 
called a "subject" 903 hereinafter. Then. - t - " stem is 

,\-:V: u-nu the X-rays and detecting the radiated X-rays, the cone beam X-ray 
CT device executes control so that the X- — ray source 901 and the plane sensor 
902 together rotate around the subject 903. Alternately, it is possible to 
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| constitute the cone beam X-ray CT device so that the X :: ray source 901 and 

the plane sensor 902 are fixed and instead the subject 903 itself rotates. 

As above, when one rotation from 0° to 360° is executed. ~ t vh *-4 the 
X- — ray source 901 and the plane sensor 902 each stands again hs >'< m- >r 
b tysrv: u>o initial positions. Alternately, when the subject 903 rotates from 0° 

to 360°, the subject 903 faces the initial direction at the initial position. 

Incidentally, at that time, the projected image initially acquired by the 
plane sensor 902 at the angle 0° has to be the same as the projected image 
acquired after the rotation of 360°. 
1 0 However, in She-actual radiography, there is no assurance that a patient 

who is b emg-the subject 903 does not move. That is, if the patient moves 
during the radiography, the projected image initially acquired by the plane 
| sensor 902 at the angle 0° is insultingly consequently d ifferent from the 
projected image acquired after the rotation of 360°. Therefore, if the 
1 5 tomographic image of the subject 903 is created by using such different 

projected images, a streak artifact appears on the created tomographic image. 

To reduce the artifact which appears due to such body movement, in 
^--conventional X-ray CT device in which the one-dimensional sensor is used, 
a method of correcting for deviation or misregistration of the subject occurring 
2 0 due to its body movement by executing interpolation within a range of certain 

angles is known (for example, see Japanese Patent Application Laid-Open No. 
H06-1 14052). 

Subsequently, this method will be explained with reference to Figs. 
| 10A, 10B and IOC. Here, Figs. 10A to 10C show fe--ajnethod -i\ nKo*: 
2 5 the body movement is corrected by using a sinogram in which the channels of 

the detector are plotted on the axis of abscissa and rotation angles are plotted on 
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the axis of ordinate. Here, it should be noted that each of the curved lines 
1001 to 1003 on this sinogram indicates the trajectory of a certain point within 
the subject 903. 

If the subject 903 does not move during the radiography (that is, there 
is no body movement of the subject 903), as shown in Fig. 10A, a detector 
channel position A of the point at the time of a start of scan (0°) conforms to a 
detector channel position B on the point at the time of an end of the scan (360°). 
On the other hand, if the body movement of the subject 903 occurs, as shown in 
Fig. 10B, the detector channel position of the point at the time of the end of the 
scan (360°) which should be essentially the position B shifts to a position B'. 

To correct for such deviation or misregistration of the subject 903 
occurring due to its body movement, as shown in Fig. 10C, the detector channel 
position A of the point at the time of the start of the scan (0°) is shifted to a 
position A and the detector channel position B' on the point at the time of the 
end of the scan (360°) is likewise shifted to a position B" both by an amount 
half as much as a deviation amount between the detector channel position A of 
the point at the time of the start of scan (0°) and the detector channel position 
B' on the point at the time of the end of the scan (360°), so as to conform the 
position A and the position B" to each other. Then, an interpolation process is 
executed te-on the data between the position A and the position B" with respect 
to all the channels to acquire the data string as- indicated by the dotted curved 
line 1003 shown in Fig. IOC. 

By using the above method, in ^-conventional X-ray CT device in 
which the one-dimensional sensor is used, it is possible to reduce the artifact 
which appears due to the body movement of the subject 903. However, in the 
cone beam X-ray CT device in which the plane sensor 902 is used, if the above 
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method is applied with respect to each of the horizontal lines constituting the 
plane sensor 902, there is a problem in that consistency between the adjacent 
horizontal lines cannot be sufficiently attained. \ m the body 

movement of the subject cannot be accurately corrected rorfe^tft^v. 
5 Moreover, if it is judged by using the magnitude of deviation or 

misregistration of each horizontal line whether or not the body movement of the 
subject has occurred*, there is a problem that accurate judgment for the body 
movement cannot be attained. . . •.. ho- !••• there is a s<>8uota h\u 

that ike-body movement >, nui; !>o needlessly corrected ».<s> >\«h 

10 

DISCLOSURE OF THE INVENTION 

The present invention has been made in consideration of the above 
problems, and the object thereof is to perf'i rm t high-accuracfely correct:- m j ir 
a body movement as much as possible, for reducing an artifact which appears 
15 on a tomographic image. 

A radiation imaging device according to the present invention is 
characterized by comprising: 

a geometric transformation parameter solving unit adapted to acquire, 
from among plural projected images of which projected angles of a radiation 
2 0 are different from others, geometric transformation parameters between at least 

the two projected images of which the projected angles of the radiation overlap 
each other; 

a changing unit adapted to gradually change the geometric 
transformation parameters acquired by the geometric transformation parameter 
2 5 solving unit, within a predetermined range of the projected angles of the 

radiation; and 
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a body movement correction unit adapted to execute a correction of a 
body movement by executing geometric transformation to the plural projected 
images of which the projected angles of the radiation are different, by using the 
respective changed geometric transformation parameters. 

Moreover, an image processing method according to the present 
invention is characterized by comprising: 

a geometric transformation parameter solving step of acquiring, from 
among plural projected images of which projected angles of a radiation are 
different from others, geometric transformation parameters between at least the 
two projected images of which the projected angles of the radiation overlap 
each other; 

a changing step of gradually changing the geometric transformation 
parameters acquired in the geometric transformation parameter solving step, 
within a predetermined range of the projected angles of the radiation; and 

a body movement correction step of executing a correction of a body 
movement by executing geometric transformation to the plural projected 
images of which the projected angles of the radiation are different, by using the 
respective changed geometric transformation parameters. 

BRIEF DESCRIPTION OF THE DRAWINGS 

Fig. 1 is a block diagram showing one example of the whole 
constitution of an X-ray CT system according to the pi-'ietuM embodiment of 
the present invention; 

Fig. 2 is a flow chart for explaining one example of the operation to be 
executed by the X-ray CT system in case of radiographing a subject and thus 
acquiring a tomographic image, according to the embodiment of R| \ to 
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f*^* : *:f-4*!V«ft$-i*m-; 

Fig. 3 is a flow chart for explaining one example of the operation to be 
executed in a pre-process and a reconstruction process of the X-ray CT system, 
according to the embodiment of ! ^ w+* \ 

Fig. 4 is a conceptual diagram for explaining the corresponding points 
between two projected images of which the projected angles overlap each other, 
according to nui <h- embodiment --*--i*K-i^iv^H?*--i-B-vV* ! *H->^; 

Fig. 5 is a conceptual diagram for explaining the fixed points on one of 
the two projected images of which the projected angles overlap each other, 
according to the embodiment of ' th e ' pr e sent invention ; 

Figs. 6A and 6B are conceptual diagrams for explaining how to acquire 
the corresponding points by means of matching, according to the foregoing 
embodiment-0f--th«--pi : e^iU--kvve«tk>ft; 

Fig. 7 is a conceptual diagram for explaining geometric transformation 
parameters for correcting body movement, according to the foregoing 
embodiment of the present invention; 

Fig. 8 is a conceptual diagram for explaining the coordinates which 
show the status before the geometric transformation is executed and the 
coordinates which show the status after the geometric transformation v*> has % - i 
executed; 

Fig. 9 is a conceptual diagram for explaining a conventional cone beam 
X-ray CT device; and 

Figs. 10A, 10B and IOC are conceptual diagrams for explaining 
conventional geometric transformation parameters for correcting for body 
movement. 
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BEST MODE FOR CARRYING OUT THE INVENTION 

Hereinafter, the ^ ( embodiment of the present invention will be 
explained with reference to the attached drawings. 

Fig. 1 is a block diagram showing one example of the whole 
5 | constitution of an X-ray CT system according to the preferred e mbodiment of 
the present invention. 
I In Fig. 1, m X-ray- ajj is generated from an X-ray source (or X-ray 

generator) 103 which is controlled by an X-ray generator control unit 104, the 
generated X rayN no ^ transmitted through a patient 102 v*v>- K >\K> is 
10 treated as a subject, and the transmitted X-rays, are te-then detected by an X- — 

ray detector 101. Incidentally, in the present embodiment, the case where *fee~ 
X-ray s are -fe-used will be explained by way of example. However, it is 
needless to say that the radiation to be processed in the present embodiment is 
not limited to lfee--X-rays-,-i that is, electromagnetic waves and such as sm- a-ray.. 
1 5 beams , a-P-ray beams , *ey-rays and the like may be used. 

The X-rays detected by the X-ray detector 101 i-s ti then input as a 
projected image to an image input unit 105. Incidentally, the X-ray source 
103 and the X-ray detector 101 such as a plane sensor or the like collect the 
projected image with respect to each of predetermined rotation angles as 
2 0 together rotating around the patient 102. However, it is alternatively possible 

to rotate the patient 102, as the subject, who is fixed to a rotary table, as 
maintaining the positional relation of the patient 1 02 with respect to the X-ray 
source 103 and the X-ray detector 101. 

Here, an image processing unit 107 corrects the X-ray detector 101, 
2 5 | executes a pre-process including logarithmic transformation to-on the input 
projected image of each of the rotation angles, and further executes image 
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processes such as a reconstruction process and the like to on the input projected 
image. Thus, a group of tomographic images (also er called a ''tomographic 
image group") is created through the above processes. Aftef tjhen, the 
created tomographic image group can be displayed on a diagnostic monitor 109, 
and also stored in an image storage unit 108. In addition, it is possible to 
output the created tomographic image group to a printer 112, a diagnostic work 
station 113 or an image database 114 through a network 111. 

Incidentally, a tester or a diagnostician handles an operation unit 110 to 
execute various operations such as a display window operation on the 
diagnostic monitor 109 or the like, a switching display operation of the 
tomographic image in the body axis of the patient 102, a cross section 
transformation operation of the tomographic image, a three-dimensional surface 
display operation, and the like. In the present embodiment, the above 
operations are totally managed by an X-ray radiography system control unit 
106 which consists of a microcomputer and the like. 

Fig. 2 is a flow chart for explaining one example of the operation to be 
executed by the X-ray CT system in case of radiographing the subject (CT 
radiography) and thus acquiring the tomographic image. 

Initially, in a step S201, radiography conditions such as an X-ray tube 
voltage, an X-ray tube current, an exposure time and the like, patient 
information such as aihe patient's name. age gender and 

the like, and test information such as a test ID (identification) and the like are 
set through the operation by the tester on the operation unit 110. 

Next, in a~step S202, the X-ray generator control unit 104 controls the 
X-ray source 103 in accordance with the radiography conditions set in -Mep 
S201. Thus, the X-ray is generated from the X-ray source 103 to radiograph 
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the patient 102, and the projected image acquired by the radiography is input to 
the image input unit 105. 

Subsequently, in «~step S203, the image processing unit 107 executes 
the pre-process te- on the projected image acquired in the radiography Then, 
5 in e-step S204, the image processing unit 107 further executes the 

reconstruction process to the projected image subjected to the pre-process in fee 
step S203. By doing so, the tomographic image of the patient 102 is created. 

After then, in a step S205, the tomographic image created in tfes-step 
S204 is displayed on the diagnostic monitor 109. Thus, the tester confirms the 
10 tomographic image displayed on the diagnostic monitor 109, and handles the 

| operation unit 110 in accordance with the confirmed result. Subsequently, in * 
step S206, the X-ray radiography system control unit 106 judges whether or not 
| it is indicated that the X-ray radiography has succeededs, based on the content 
of the handling by the tester. Then, when the X-ray radiography system 
15 control unit 106 judges that it is indicated that the X-ray radiography has d ees- 

not succeeded, the flow returns to fee-step S202 to a-; «*-.iyj a? radiography of 
the patient 102. That is, the processes in fee-steps S202 to S206 are repeated 
until the X-ray radiography system control unit 106 judges that it is indicated 
| that the X-ray radiography has succeededs. On the contrary, when the X-ray 
2 0 radiography system control unit 106 judges that it is indicated that the X-ray 

| radiography succeeded*, the flow advances to a-step S207. In fee-step S207, 
the X-ray radiography system control unit 106 transfers the tomographic image 
created in fee-step S205 to the image storage unit 108, the image database 114 
and the like. Subsequently, in e-step S208, the X-ray radiography system 
2 5 control unit 106 judges whether or not next radiography is indicated, on the 

basis of the handling by the tester on the operation unit 110. Then, when the 
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X-ray radiography system control unit 106 judges that the next radiography is 
indicated, the flow returns to fee- step S201 to execute the next radiography. 
On the contrary, when the X-ray radiography system control unit 106 judges 
that the next radiography is not indicated, the operation to be executed by the 
X-ray CT system ends. 

Fig. 3 is a flow chart for explaining one example of the operation to be 
executed in the pre-process and the reconstruction process of the X-ray CT 
system, according to the present embodiment. 

In Fig. 3, at the first place, the coordinates of the corresponding points 
between the two projected images of which the projected angles overlap each 
other are acquired (step S301). Here, in the following, a concrete example of 
the process to be executed in {-in^-step S301 will be explained. 

More specifically, for example, if the set of the X-ray source 103 and 
the X-ray detector 101 scans the patient 102 as just making a round (360°) of 
the patient 102, the projected angles which overlap each other may be set as 0° 
and 360° respectively. \ ,M v?*4^ , if the set of the X-ray source 103 and the 
X-ray detector 101 executes the scan of 360° or more ((360 + a)°), the projected 
angles which overlap each other may be set as 9° (< a) and (360 + 0)°, 
respectively. 

M Fig. 4 is a conceptual diagram for explaining the corresponding 
points between the two projected images of which the projected angles overlap 
each other. That is, as shown in Fig. 4, the coordinates of the corresponding 
points (the small black square points shown in Fig. 4) between a projected 
image 401 (for example, an image at scan angle 0°) and a projected image 402 
(for example, an image at scan angle 360°) of which the respective projected 
angles overlap each other are acquired. Here, it should be noted that the 
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number of corresponding points has to be equal to or larger than the number of 
geometric transformation parameters, and accuracy of the geometric 
| transformation parameters improves in proportion to an.increase in the number 
of corresponding points. 
5 However, if the number of corresponding points increases, athe 

necessary calculation time becomes long > \ ^sv,^ 

Therefore, in consideration of such a fact, it has only to determine the number 
of corresponding points substantially in accordance with the system itself. In 
any case, the numbei oj sets of the coordinates of the respective corresponding 

10 points of the projected images 401 and 402 : - are acquired .is ^\ -s 

equal to the number of corresponding points. To achieve this, first, plural 
fixed points are set on one (e.g., projected image 401) of the two projected 
images. Here, for example, it is preferable to use the most characteristic 
points as the fixed points. For example, with respect to the image of a chest, 

15 the end points of a rib, the branch points of lung blood vessels, and the like are 

suitable for the most characteristic points. Therefore, the end points of the rib, 
the branch points of the lung blood vessels or the like are set as the fixed points. 
In contrast, even if the characteristic points are not used as the fixed points, for 
example, as shown in Fig. 5, it has only to divide the entirety of the projected 

2 0 image 401 into numbers of small rectangles, and then set the respective central 

points of the divided small rectangles within the range of the subject as the 
fixed points. 

Thus, if the plural fixed points are set with respect to one of the two 
projected images, the coordinates, on the other (e.g., projected image 402) of 
25 the two projected images, of the points respectively corresponding to these 

fixed points are acquired. To achieve this, it has only to use a matching 
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method which is generally used in stereo analysis and the like. Hereinafter, an 
example of how to acquire the sets of the corresponding points by using the 
matching method will be explained with reference to Figs. 6A and 6B. 

As shown in Fig. 6A, a template image 603 of a predetermined size is 
cut out on the basis of the fixed point on the projected image 401. Then, as 
shown in Fig. 6B, a search area 604 of a predetermined size larger than the size 
of the template image 603 is set on the projected image 402 of which the 
projected angle overlaps the projected angle of the projected image 401. Here, 
it should be noted that the coordinates of the central point of the search area 
604 are the same as the coordinates of the fixed point on the projected image 



Then, the position of the central point is shifted sequentially from the 
upper left to the lower right within the search area 604. At the same time, a 
reference image of the size same as that of the template image 603 is cut out 
from the projected image 402. Subsequently, an estimated value est 
concerning the matching between the cut-out reference image and the template 
image 603 is acquired. 

Here, if it is assumed that the coordinates of the reference image are 
given as r(x, y) and the coordinates of the template image 603 are given as t(x, 
y), the estimated value est of the matching can be acquired by the following 
equation 




est - 
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Here, it should be noted that, in the above equation (1), the symbols n 
and m respectively denote pixel sizes of the template image 603, and the 
horizontal lines (-) added above the characters fr and tj indicate the i^eettw 
average values. . ; ; • -i . 

5 Incidentally, when it is requested to execute the calculation at high 

speed, the estimated value est of the matching may be acquired simply by the 
| following equation (2) :-- 

<*'=S2>(*,/)-f(*,/)| 

10 

That is, the position of the reference image at which the estimated 
value est of the matching becomes minimum is equivalent to the coordinates of 
the corresponding points. 

In any case, when the coordinates of the corresponding points between 
15 the projected images 401 and 402 of which the projected angles overlap each 

| other are acquired in fee-step S301, the geometric transformation parameter is 

acquired from the set of the coordinates of the corresponding points (step S302). 
| For example, an affine transformation may be used in such geometric 
transformation. 

2 0 It should be noted that the affine transformation can be represented by 

| the following equation (3)< 




25 

Then, the set of the corresponding points is defined as the set of the 
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point of the coordinates (x„, y n ) and the point of the coordinates (X„, Y„) (n is 
1, N). Subsequently, geometric transformation parameters a 0 i, a 0 2, an, a J2 , 
a 2 i, a 2 2 for the affine transformation by which a least squares error err between 
the coordinates (x n \ y n ') acquired by transforming the coordinates (x„, y„) with 
use of the above equation (3) and the coordinates (X n , Y„) is minimized are 
acquired. Incidentally, the least squares error err is represented by the 
following equation (4):.- 



Here, the coordinates x n ' and y n ' are the functions which use the 
geometric transformation parameters a 0 i, a 0 2, an, a ]2 , a 2 i, a 22 for the affine 
transformation as variables. Therefore, also the least squares error err by the 
equation (4) is the function which uses the geometric transformation parameters 
a 0 i, a 02 , an, ai 2 , a 2 i, a 22 for the affine transformation as variables. 

For this reason, if the least squares error err is subjected to partial 
differentiation based on each of the variables of the geometric transformation 
parameters aoi, ao 2 , an, aj 2 , a 2 i, a 22 for the affine transformation to set "0" and 
then the simultaneous equations are solved, it is possible to acquire the 
geometric transformation parameters a 0 i, a 02 , an, ai 2 , a 2 i, a 22 for the affine 
transformation by which the least squares error err becomes minimum. 

In this way, when the affine transformation is executed, the geometric 
transformation parameters a 0 i, ao 2 , an, ai 2 , a 2 i, a 22 for the affine transformation 
are acquired. However, when the affine transformation is executed, the 
geometric transformation in which also enlargement of the image and shearing 
of the image are admitted. Consequently, in case of executing the geometric 




15 Substitute Specification (marked) 

A.N. 10/599,028 
Atty. Docket No. 03500.109704. 

transformation only for rotation of the image and movement of the image 
without permitting the enlargement and shearing of the image, 

Helmert transformation is used. Here, it should be noted that the Helmert 
transformation can be represented by the following equation (5)~; 

::km:} - 

That is, the geometric transformation parameters in case of the Helmert 
transformation are three, i.e., 9, a 01 , a 02 . Therefore, these parameters can be 
acquired in the same manner as above. 

In addition to the affine transformation and the Helmert transformation, 
even if a secondary projective transformation, high-order polynomial 
transformation or the like is used, it is possible to acquire the geometric 
transformation parameters as well as the case where the affine transformation, 
the Helmert transformation or the like is used. 

In any case, when the geometric transformation parameter is acquired 
in the-step S302, it is then judged whether or not an estimated amount of the 
acquired geometric transformation parameter is larger than a predetermined 
amount (step S303). 

Here, it should be noted that the estimated amount of the geometric 
transformation parameter is the amount for estimating the magnitude of the 
body movement while the scan is being executed. More specifically, the 
estimated amount of the geometric transformation parameter is used to estimate 
how far from the geometric transformation parameter in case of no body 
movement. 

Incidentally, for example, it is assumed that the geometric 
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transformation parameters are set as Xi, X 2 , —, K, —, hi, and the geometric 
transformation parameters in a.case of no body movement are set as V, 
Xn, —, ^n'- On such a premise, an estimated amount est2 of the geometric 
transformation parameter can be acquired by the following equation (6);- 

Here, it should be noted that the right side of the above equation (6) 
may not be necessarily the square root. That is, instead of the above equation 
(6), it is possible to more simply acquire the estimated amount est2 of the 
geometric transformation parameter by the following equation (7);- 

est2 = f}A'„-X„\-(7) 

M = I 

By the way, if the affine transformation is adopted as the geometric 
transformation, the following equations (8) are satisfied.^ 
X\ = aoi, A,2 = ao2, ^3 = an> ^4 = &i2> ^.5 = a 2i> 
= a 22 

v = o, 12 = o, v = 1, w = 0, v = 0, v = 1 

... (8) 

Moreover, if the Helmert transformation is adopted as the geometric 
transformation, the following equations (9) are satisfied:- 
A4 = 9, 7yi = aoi, A,3 = ao2 

V = 0,V = 0,V = 0 -(9) 

In any case, when it is judged in fee -step S303 that the estimated 
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amount of the geometric transformation parameter acquired as above is equal to 
or smaller than a predetermined amount, it is possible to judge that there is no 
body movement of the level for which the correction is necessary. Therefore, 
in this case, the flow advances to a-step S306 without executing the correction 
5 t»£-for the body movement so as to reconstruct the projected image and thus 
create the tomographic image. 
; On the other hand, when it is judged in the -step S303 that the estimated 

amount of the geometric transformation parameter acquired as above is larger 
| than the predetermined amount, the correction ef-for the body movement is 
1 0 executed. 

More specifically, to execute the correction e£-for the body movement, 
the geometric transformation parameter for geometric correction is first 
determined so that the projected images of which the projected angles overlap 
each other conform to each other and gradually change within a predetermined 
15 angle (step S304). 

■ In the correction ef -for the body movement, only half of the correction 

amount is gradually subjected to the correction according to the angles within 
the range, e.g., from 0° to 180°. Subsequently, within the range from 180° to 
360°, the correction is gradually executed according to the angles in the 
2 0 direction opposite to the direction of the correction within the rage from 0° to 

180°, whereby it is possible to achieve the smooth correction of the body 
movement. 

With respect to the geometric transformation parameters for executing 
the above correction of the body movement (geometric correction), only half of 
25 | the geometric transformation parameters acquired in fee-step S302 may be 
determined so as to gradually change within the range, e.g., from 0° to 180°. 
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Besides, only half of the minus geometric transformation parameters acquired 
| in the-step S302 may be determined so as to gradually change within the range, 
e.g., from 180° to 360°. 

In addition, the body movement is gradually corrected according to the 
5 angles within the range from 0° to 90°, by half of the correction amount of the 

body movement. Within the range from 90° to 270°, the body movement is 
not corrected for. Then, within the range from 270° to 360°, the body 
movement is gradually corrected for according to the angles, by half of the 
correction amount of the body movement, in the direction opposite to the 
10 direction of the correction within the rage from 0° to 90°. 

More specifically, the meaning that the projected image gradually 
changes within the predetermined angle and the body movement is gradually 
| corrected for according to the angles is as follows. That is, the geometric 
transformation parameter is divided by the predetermined angle to acquire a 
1 5 change amount of the geometric transformation parameter with respect to each 
unit angle, and the change amounts of the geometric transformation parameters 
of the respective unit angles are gradually added to others. 

Here, one example of the method of determining the geometric 
transformation parameters for executing the correction of the body movement 
2 0 (geometric correction) will be explained with reference to Fig. 7. 

That is, Fig. 7 shows the projected angles in the direction along the 
axis of abscissas and the geometric transformation parameters in the direction 
along the axis of ordinate, with respect to each of the respective geometric 
transformation parameters. In Fig. 7, the angles from 0° to 360° are given as 
25 the projected angles. However, when the scan of an angle exceeding 360° 

((360 + a)°) is executed, the projected angle may be set from 9° (< a) to (360 + 
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9)°. 

Here, it takes the case of n-th parameter as an example. The 
geometric transformation is executed by half of the body movement with 
respect to the projected image of 0°. whereby a nd so its parameter n„(0) is 
5 given by the following equation (10):.- 

Mn(0) = V + (k-V)/2 ...(10) 

Then, for example, it is assumed that the above predetermined angle is 
of. In such a case, since the geometric transformation is not executed within 
the range from co° to (360 - co)°, the parameter of the projected angle within this 
1 0 range is assumed as A. n ' for the sake of convenience. Here, within the range 

from 0° to co°, it is controlled to gradually change the parameter from n„(0) to 
X, n '. That is, the parameter u„(P) in case of P° (< co) is given by the following 
| equation (11 ) : ;. 

^n(P) = V+{an-V)/2}xP/(0 ...(11) 

1 5 Similarly, the geometric transformation is executed by half of the body 

movement with respect to the projected image of 360°, wheffrbv ^ir..:. its 
parameter |a n (360) is given by the following equation (12).:.-- 
Hn(360) = K'-(K- V)/2 ...(12) 

Then, within the range from (360 - co)° to 360°, it is controlled to 
2 0 gradually change the parameter from A, n ' to u. n (360). That is, the parameter 

|i n (y) of the range from (360 - co)° to 360° is given by the following equation 
| (13), 

Hn(y) = K' +{(K- K')/2] x {y - (360 - Q))}/(0 

... (13) 

2 5 These parameters are respectively used to execute the geometric 

correction by half with respect to the side of 0° and also to execute the 
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geometric correction by half with respect to the side of 360°. However, the 
parameter may be determined so that the whole transformation is executed 
within the predetermined angle with respect to either the side of 0° or the side 
of 360°. 

5 Then, when the geometric transformation parameters for correcting for. 

the body movement are acquired, the geometric transformation is executed te- 
em the projected image of each of the projected angles by using the acquired 
geometric transformation parameter (step S305). However, it should be noted 
that the geometric transformation is not executed within the range from co° to 

10 (360-co)°. 

; Now Subseque mly. one example of the geometric transformation will 

be explained with reference to Fig. 8. In Fig. 8, it is assumed that, through the 
geometric transformation, a point (X, y) indicated by numeral 803 on a 
projected image 801 which is acquired before the geometric transformation is 

15 | executed (hereinafter, called a. "before-geometric transformation imaged 801) is 
transformed into a point (X\ y') indicated by numeral 804 on a projected image 
802 which is acquired after the geometric transformation is executed 
| (hereinafter, called an "after-geometric transformation image" 802). 

When an image is created through the geometric transformation, the 

2 0 process opposite to the above is used. That is, the point (X, y) on the before- 

geometric transformation image 801 to be transformed into the point (X', y') on 
the after-geometric transformation image 802 is acquired. Here, it should be 
noted that the point (X', y') indicates the coordinates of the integer value 
corresponding to each pixel on the after-geometric transformation image 802. 

2 5 Then, by using the inverse transformation of the geometric transformation with 

respect to the coordinates, the point (X, y) corresponding to the point (X', y') 
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can be acquired. However, since the point (X, y) does not have the 
coordinates of an integer value, interpolation is executed by using the pixel 
values of the four points closest to the point (X, y), and the interpolated value is 
set as the pixel value of the point (X', y'). This process is executed te-on all 
the pixels, whereby the image can be created through the geometric 
transformation. 

Incidentally, when the image is created through the geometric 
transformation as above, the geometric transformation is not directly used, but 
the inverse transformation of the geometric transformation is used. For this 
reason, the geometric transformation parameter may be set in advance as the 
parameter for the inverse transformation of the geometric transformation in #*e- 
step S304. 

In any case, when the correction ef-for the body movement through the 
geometric transformation ends, the reconstruction is executed by using the 
corrected projected images, and the tomographic image is thus created (step 
S306). 

As described above, in the present embodiment, the coordinates of the 
respective corresponding points between the projected images 401 and 402 of 
which the respective projected angles (for example, 0° and 360°) overlap each 
other are first acquired. Then, the geometric transformation parameters |i n (0) 
and |i n (360) are acquired through the affine transformation or the like by using 
the sets of the coordinates of the corresponding points. Subsequently, when 
the estimated amount of the geometric transformation parameter is larger than 
the predetermined amount, the correction ef-for the body movement is executed, 
and the tomographic image is then created by using the projected image 802 
which has been subjected to the correction ef-for the body movement. 
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More specifically, when the correction ef -for the body movement is 
executed, as described above, the geometric transformation parameter for 
executing the geometric correction with respect to the projected angles from co° 
to (360 - (o)° is assumed as V for the sake of convenience. Then, the 
geometric transformation parameter |i n ((3) for executing the geometric 
correction is determined so that the parameter gradually changes from ^i n (0) to 
X n ' (as indicated by the equation (11)) within the range of the projected angles 
from 0° to cd°. Moreover, the geometric transformation parameter n„(y) for 
executing the geometric correction is determined so that the parameter 
gradually changes from to |a n (360) (as indicated by the equation (13)) within 
the range of the projected angles from (360 - co)° to 360°. After then, the 
geometric transformation (i.e., body movement correction) is executed te-on the 
projected images having the projected angles from 0° to co° and te -on the 
projected images having the projected angles from (360 - co)° to 360°. 

As described above, according to the present embodiment, the body 
movement (that is, the movement of the subject) is not grasped all -\ j with 
respect to each horizontal (landscape) line, but the body movement is gra ^pe**- 
as the image. Thus, the body movement can be accurately uU ;u:ivJ mj 
<>v ■ >;:• -b'.i'i-U '!>:•! v-i-/e-d, whereby it is possible to achieve high-precision 
correction eft-for the body movement, and it is also possible to achieve the- 
smooth correction ef-for the body movement. In consequence, it is possible to 
create f^e :; tomographic image on which an artifact has been reduced by 
correcting for the body movement. 
(Other Embodiments Of The Present Invention) 

Incidentally, the present invention includes in its category a case where 
; : program codes of software for realizing the functions of the above 
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| embodiment is m supplied to a computer provided in an apparatus or a system 
connected to various devices so as to actually operate these devices for 
realizing the functions of the above embodiment, and thus the computer (CPU 
or MPU) in the system or the apparatus reads and executes the supplied and 
5 stored program codes and operates these devices. 

In this case, the program codes themselves of the software realize the 
functions of the above embodiment. Therefore, the program codes themselves, 
and a means? such as a recording medium storing these program codes, for 
supplying the program codes to the computer constitute the present invention. 
10 As the recording medium for storing these program codes, for example, a 

flexible disk, a hard disk, an optical disk, a magnetooptical disk, a CD-ROM, a 
magnetic tape, a nonvolatile memory card, a ROM or the like can be used. 

Moreover, it is needless to say that the present invention includes not 
only a case where the functions of the above embodiment are realized by 
1 5 | executing the supplied program codes with the y sinjij. computer, but also a 

case where the program codes cooperate with an OS (operating system) running 
| on the computer or another application software to realize the functions of the 
above embodiment. 

Furthermore, it is needless to say that the present invention also 
2 0 includes a case where, after the supplied program codes are written into a 

function expansion board inserted in the computer or a memory of a function 
expansion unit connected to the computer, the CPU or the like provided in the 
function expansion board or the function expansion unit executes a part or all of 
the actual processes on the basis of the instructions of the program codes, and 
25 thus the functions of the above embodiment are realized by such the processes. 

According to the present embodiment, She-body movement of the 
c-n^i vt can be determined from g rasped as t he image. Thus, the body 
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| movement can be accurately Mom !tic« I nso «n>ic«l<f>- x In 

addition, the geometric transformation parameters are smoothly changed 
between the projected angles, whereby it is possible to achieve high-precision 
correction ^1 u N r the body movement as much as possible. In consequence, it 
5 is possible to create fe*-a tomographic image on which an artifact has been 

reduced as much as possible. 

According to another aspect of the present invention, the magnitude of 
the body movement is estimated by using the acquired geometric 
| transformation parameters, and the correction e£ -for the body movement is 
1 0 executed only in the case where the magnitude of the body movement is equal 

to or larger than the predetermined magnitude. Thus, useless correction of the 
body movement can be avoided, whereby it is possible to shorten the time 
| necessary for executing the correction eg-forthe body movement. 

1 5 This application claims priority from Japanese Patent Application No. 2004- 
141490 filed May 11, 2004, which is hereby incorporated by reference herein. 
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ABSTRACT 

To ivrtorm high-accurate -\ .-k 0 > ^..^ as possible a 
correction ef for a body movement for reducing an artifact occurring in a.,case 
of creating a tomographic image from an X-ray projected image, the 
5 | coordinates of the respective corresponding points are acquired as_between the 
projected images of which the projected angles overlap each other (e.g., 0° and 
| 360°), and #*e-a geometric transformation parameter is acquired through affine 
transformation or the like by using the set of the acquired coordinates. When 
an estimated amount of the acquired parameter is equal to or larger than a 
1 0 predetermined amount, the geometric transformation parameter for geometric 

correction is determined by using the acquired geometric transformation 
parameter, the geometric transformation (correction of body movement) is 
executed by using the determined geometric transformation parameter, and the 
tomographic image is created by using the corrected projected images. 
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